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LONG-TERM  GOAL 

To  develop  methods  to  utilize  synthetic-aperture  radar  (SAR)  to  characterize  wave  coherence  on  scales 
relevant  to  the  design  of  mobile  off-shore  bases  (MOB). 

OBJECTIVES 

There  are  two  basic  objectives  to  this  program:  (1)  to  develop  the  appropriate  measures  of  wave 
coherence  and  methods  of  application  for  SAR  image  data,  and  (2)  to  apply  the  methods  to  available 
SAR  data  sets. 

APPROACH 

Methodology  for  determining  the  crest-length  distribution  using  SAR  image  data,  including  error 
metrics,  were  developed  and  then  validated  using  available  ground  truth. 

WORK  COMPLETED 

Methods  for  characterizing  crest-length  distributions  have  been  developed  and  applied  to  both 
simulated  and  actual  SAR  data.  The  impact  of  SAR  imaging  effects  on  the  results  for  swell-like 
conditions  has  been  assessed  and  guidelines  have  been  developed  to  mitigate  that  impact.  For  more 
complicated  spectra,  a  variational  method  for  estimating  the  wave  spectrum  from  the  SAR-image 
spectrum  has  been  developed. 

RESULTS 

The  results  of  this  study  are  in  two  parts.  First  methods  for  assessing  wave-coherence-related 
quantities  directly  from  SAR-image  data  for  simple  swell-like  spectra  were  developed.  For  more 
complicated  situations,  a  variational  inversion  algorithm  was  developed  to  determine  the  wave 
spectrum  from  the  SAR-image  data.  The  estimated  wave  spectrum  is  then  used  to  characterize  the 
wave  coherence.  The  key  elements  of  these  two  approaches  are  presented  below. 

The  distribution  of  wave-crest  lengths  on  the  ocean  surface  can  be  shown  to  primarily  depend  on  the 
directional  spread  of  the  wave  spectrum.  It  is  well  known,  however,  that  SAR  imaging  effects  can  alter 
the  apparent  directional  spread  of  the  wave  spectrum.  How  the  SAR  imaging  effects  manifest 
themselves  depends  on  the  wavelength  of  the  waves,  the  significant  wave  height,  and  the  SAR  look 
direction  relative  to  the  wave  propagation  direction.  To  investigate  the  effects  of  SAR  imaging  on  the 
apparent  crest-length  distribution  derived  from  a  SAR  image  for  swell-like  conditions  (i.e.  for  a 
unimodal  wave  spectrum,  relatively  Gaussian  in  shape),  a  series  of  forward  predictions  of  the  SAR 
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Figure  1:  Apparent  directional  spread  Atj) 
from  simulated  SAR  image  normalized  with 
actual  directional  spread  A<p  =  0.105  for 
various  SAR  look  directions  relative  to  the 
dominant  wave  direction  for  H,  =  2.0  m,  5  m 
and  10.0  m  for  X  =  500  m. 


Actual  Wave  Direction  (degrees  from  SAR  look  dir) 

Figure  2:  Apparent  wave  direction  for 
A<p  Actual  =  0.105  VS.  actual  wave  direction  with 
=  2.0  m,  5  m  and  10.0  m,  and  X  =  500  m. 


image  spectrum  for  known  wave  spectra  were  carried  out  using  the  Hasselmann  formulation  [1,2]. 
Simulations  were  made  over  a  range  of  significant  wave  heights  (2-10  m),  dominant  wave  frequencies 
(0.06-0.09  Hz),  and  SAR  look  directions  relative  to  the  wave  propagation  direction.  The  SAR 
parameters  for  the  simulations  were  those  typical  of  satellite  SAR  sensors  (specifically  the  European 
Space  Agency’s  ERS  platform). 

Eigure  1  shows  the  apparent  directional  spread  of  the  SAR  spectrum,  divided  by  the  directional  spread 
of  the  actual  wave  spectrum,  as  a  function  of  the  SAR  look  direction  relative  to  the  dominant  wave 
direction  for  three  significant  wave  heights.  These  results  are  for  a  narrow  wave  spectrum  {A(jt=0.\05 
rad)  and  a  dominant  wavelength  of  500  m.  Eor  all  wave  heights,  the  error  increases  with  look  direction 
and  becomes  constant  at  an  apparent  directional  spread  which  is  about  30-40%  too  low  near  a  relative 
look  angle  of  20°.  The  apparent  narrowing  of  the  wave  spectrum  is  due  to  the  narrow  azimuth  pass 
band  associated  with  the  SAR  imaging  process,  and  the  consequent  distortion  of  the  wave  spectrum. 
The  errors  for  this  case  are  conservative  in  the  sense  that  the  apparent  directional  spread  would  indicate 
that  the  crests  longer  than  they  actually  are.  Eor  shorter  wavelengths,  and  larger  directional  spreading  , 
the  effects  can  be  more  pronounced. 

Erom  figure  1,  it  is  clear  that  the  amount  of  error  in  the  directional  spread  determined  from  the  SAR 
image  spectrum  depends  on  the  look  direction  of  the  SAR  relative  to  the  wave  direction.  The  apparent 
wave  direction  can  be  estimated  from  the  location  of  the  peak  in  the  SAR  spectrum.  Eigure  2  shows 
the  apparent  wave  direction  plotted  versus  SAR  look  direction  for  the  spectra  discussed  in  figure  1 . 
Erom  the  figure,  it  is  clear  that  the  wave  direction  can  be  estimated  within  a  few  percent  from  the  SAR 
spectrum  under  all  conditions.  Hence,  the  wave  direction  can  be  estimated  accurately  from  the  SAR 
spectrum  for  swell-like  conditions. 

Eor  more  complicated  spectra,  ignoring  the  SAR  imaging  effects  and  treating  the  SAR-image  spectrum 
as  if  it  were  the  wave  spectrum  can  lead  to  serious  errors.  As  a  result,  a  variational  inversion  scheme 
based  on  the  Hasselmann  &  Hasselmann  model  was  developed.  This  effort  leveraged  and  extended 


work  done  on  assimilation  of  SAR  data  into  the  SWAN  wave-speetrum  model  under  the  ONR 
Advaneed  Wave  Predietion  Program.  The  inversion  proeedure  involves  adjusting  the  wave  height 
speetrum  S(k)so  as  to  minimize  the  eost  funetion 

7  =  j[5,(£)-5.(£)fj£ 

where  5,(k)  is  the  predieted  image  speetrum  and  S.(k)  is  the  observed  image  speetrum.  The  predieted 
image  speetrum  is  ealeulated  from  the  equation 

=  jj  G{x,y,kje~‘''^’^'‘^'‘^^^dxdy 

where  G(x,y,kJ  is  related  to  the  wave  speetrum  as  deseribed  by  Hasselmann  and  Hasselmann  [1]  and 
Krogstad  [2].  This  relationship  involves  the  radar  modulation  transfer  funetion,  for  whieh  we  use  a 
modifieation  of  the  expressions  given  by  Plant  and  Zurk  [3].  The  ehange  in  the  eost  funetion  due  to  a 
small  ehange  in  the  wave  speetrum  ean  be  written  as 

SJ  =  ^  P(k)dS(k)dk 

where  P(k)  ~  [5,(k)  -  5,(k)][R^(k)  +  Rf(-k)]  andR^(k)is  the  quasi-linear  SAR  modulation  transfer 
funetion.  Thus,  the  eost  funetion  ean  be  redueed  by  ehanging  the  wave  speetrum  by  the  amount 
dS{k)  =  -ePik)  where  e  >  0.  Starting  with  an  initial  estimate  of  the  speetrum  (usually  zero),  we 
ealeulate  5,(k)  and  P{k) ,  and  ehange  the  speetrum  using  this  equation  until  J  is  minimized. 

Figure  3  shows  the  ehange  in  the  eost  funetion  as  a  funetion  of  the  iteration  number  for  an  example  ERS 
data  set  near  Duek,  NC.  Figure  4  shows  the  ERS  image  speetrum,  the  estimated  wave  height  speetrum 
after  50  iterations,  the  image  speetrum  eomputed  for  this  wave  speetrum,  and  the  wave  height 
speetrum  measured  at  the  FRF  8-m  array  and  projeeted  out  to  the  region  where  the  ERS  speetrum  was 
measured.  Figure  4(a)  represents  an  average  of  the  speetra  eomputed  from  16  eontiguous  128x128  pixel 
regions  of  the  ERS  image.  This  average  speetrum  was  eorreeted  for  the  stationary  impulse  response 
funetion  and  smoothed  using  a  Gaussian  kernel  with  a  width  of  two  wavenumber  samples.  The 
direetional  width  of  this  SAR-image  speetrum  is  about  twiee  that  of  the  FRF  speetrum  (figure  4d). 
Applying  the  same  measure  to  the  estimated  wave  height  speetrum  results  in  a  width  only  about  30% 
larger  than  the  FRF  speetrum  (as  projeeted  to  the  depth  of  the  ERS  measurement).  This  indieates  the 
utility  of  the  SAR  inversion  algorithm  in  improving  wave  eoherenee  estimates  from  SAR  image  data. 

IMPACT/APPLICATION 

The  development  of  approaehes  to  use  satellite-based  SAR  data  to  determine  wave  eoherenee  for  MOB 
design  will  allow  the  available  arehive  of  SAR  data,  whieh  spans  a  number  of  years,  to  be  used  to  obtain 
this  information.  This  may  eliminate  the  need  for  a  eostly  field-measurement  program  to  obtain  this 
data. 

TRANSITIONS 

Results  were  presented  at  two  MOB  Teehnology  Exehange  Conferenees  and  the  1999  Workshop  on 
Very  Farge  Floating  Struetures. 
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Figure  3.  Cost  function  and  significant  wave  height  for  wave  spectrum  estimated  from  ERS  08/27/97 

data  set,  plotted  versus  number  of  iterations. 
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Figure  4.  SAR  spectrum  inversion  results  for  ERS  SAR  data  set  collected  on  08/27/97  near  Duck, 
NC.  (a)  ERS  image  spectrum,  (b)  estimated  wave  height  spectrum,  (c)  image  spectrum  predicted  from 
estimated  wave  height  spectrum,  (d)  wave  height  spectrum  measured  at  FRF  8-m  array  and  projected 
to  approximately  20m  depth  where  ERS  image  spectrum  was  extracted. 


RELATED  PROJECTS 


This  project  is  related  to  other  efforts  under  the  MOB  Wave  Coherence  program.  These  other  efforts 
focus  on  developing  both  improved  measures  and  understanding  of  nonlinearity  and  wave  coherence. 
The  basis  for  the  variational  approach  for  estimation  of  wave  spectra  from  SAR-image  data  was 
initially  developed  under  the  ONR  Advanced  Wave  Prediction  Program  (Contract  No.  N00014-98-C- 
0012). 
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